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Abstract 

In the paper we offer a functional-discrete method for solving the Cauchy 
problem for the first order ordinary differential equations (ODEs). This 
method (FD-method) is in some sense similar to the Adomian Decomposition 
Method. But it is shown that for some problems FD-method is convergent 
whereas ADM is divergent. The results presented in the paper can be easily 
generalized on the case of systems of ODEs. 
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1. Introduction. 



Many scientific papers devoted to the Adomian decomposition method 
(ADM) have been published during last two decades. For the first time this 
method was proposed by the American physicist G. Adomian as a method for 
solving an operator equations (see Adomian 1984 -Adomian 1994] ) . The 
method is based on the specific analytical representation of the exact solution. 
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More precisely, the solution is represented in terms of a rapidly convergent 



infinite series with easily computable terms (see Cherruault. 1989 -El-Kalla 
] 2007 and the references therein). In spite of the considerable scientific 
activity, the necessary and sufficient conditions that provide a convergence 
of these series are unknown. Partially, this topic was discussed in Cherru- 
ault. 1989 - Hosseini, Nasabzadeh 2006 , where authors study the ques- 



tion about convergence of ADM applied to the nonlinear operator equation 
y — N (y) = f, where N (•) is a nonlinear operator in a Hilbert space H . The 

offers sufficient conditions 



Himoun, Abbaoui and Cherruault 2003 



paper 

that provide a convergence of ADM applied to the Cauchy problem on a fi- 



nite segment. In Hashim, Noorani, Ahmad, Bakar, Ismail and Zakaria 2006 



authors apply ADM (in a rather modified form) to the Lorenz system. The 



paper Inc and Cherruault 2005 is devoted to the application of ADM to the 



solution of the nonlinear Volterra-Fredholm integro-differential equation. 

proposed a modification of the ADM. 



The author of El-Kalla 2007 



Using several numerical examples, he has shown that this modified method 
converges faster than the ordinary one. Nevertheless, one crucial fact was 
overlooked in his work, as it turns out, this modified ADM coincides with 
the fixed point iteration method (c.f. Gavrilyuk, Lazurchak, Makarov and 
Sytnyk|[2CJ09] ). 



The idea of ADM is similar to the idea of a functional-discrete method, 

In this 



called FD-method, which was firstly introduced in Makarov 1991 



paper the FD-method was applied to the solution of the Sturm-Liouville 
problem and the remarkable convergence results were obtained. Then, in 



Gavrilyuk, Klymenko, Makarov and Rossokhata 2007 - Gavrilyuk et al 



2009 , FD-method was successfully applied to several operator equations 
and, in particular, to the boundary value problems. 

The essential difference between ADM and FD-method is expressed by 
the fact that the last one has a built-in adjustable parameter, by varying 
which we can provide the convergence of FD-method even, when ADM is 
found to be divergent. 

Let us outline the general idea of FD-method by applying it to the fol- 
lowing Cauchy problem: 



L r (u (x)) — N (x,u (x)) u(x) = (ft (x) , x G [x , +oo) 



U(x ) = lt , j-U (x ) = U ( o\ 



d r_1 / \ 



u 



(r-1) 



(1) 



2 



where L r (•) 



r-l 



dx r 



E Pk (x) 



k=l 



dx( r ~ k ) 



is a linear differential operator of 



the r-th order, r e N,p k (x) G C ([x , +00)) ; A" (x,u): [x , +00) x R — >• R 
is a given nonlinear function, 4>(x) G C ([xq, +00)) . From now on we make 
the assumption that N (x, u) is continuous in x on [x , +00) and infinitely 
differentiable with respect to u on R, i.e., iV (x, u) G C°'£° ([x , +00) x M.) . 
The main idea of FD-method is to find an approximation to the exact so- 

m 

lution u (x) of problem ( 1 ) in the form of the finite subsum u (x) = ^ u^' (x) 

i=0 

of the infinite series representation 



U (x) = ^ « W (x) ■ 



(2) 



i=0 



(3) 



To define each term of series (J2J) we need to introduce a grid: 

u = {xq < Xi < x 2 < ■ ■ ■ , x n —> +00, n — > +00} , 

h = sup {hi = Xi — Xi-i} . 

We define the functions (x) G C r_1 ([x , +00)) to be the solutions of the 
following linear Cauchy problems 

L r (x)) - AT (x, (xi_i)) w (0) (x) = (x) , 



u(°) (io) = <',±«(io)=< , l 



,(0) d 



(1) 



(x ) = U, 



(r-l) 



(4) 



x G [xj_i, Xj] , i 1,2,..., 
L r (x)) - AT (x, (x) = 

iV' (x, (x^)) (x) uCj+i) (^i) + (x) 



u 



(i+i) 



(x ) = 0, £uU+V (x ) = 0, 



(i+i) 



(5) 



x ) = 0, 



with the matching conditions 



1,2,... 



(*) 



lun_£u^ (x) 



lim 



(x) = 0, 



x= Xi x^Xi+0 " x x-^Xi-0 

fc = 0,l,...,r-l; z = l,2,...;j = 0,l,... 



(6) 



3 



where 

( x ) = £ Aj+l . p (N (x, (•)) ; t*<°> (x 4 _0 , . . . , (^0) u« (x) + 

p=i 

+ E (■)) ; « (0) (*) , (*),..., (*)) - 

p=0 

- Vp (•)) ; « (0) (^-1) , (^-1) , • • • , fo-i))] « (P) (a?) + 

(JV (x, (•)) ; (x^) , . . . , u (j) ( Xi „ x ) , 0) U W (x) , j = 0,1, — 

And 

A k (N(x, (■)); u^,u^,...,u^) 



1 d k 



i=0 



(7) 
(8) 



t=o 



are well-known Adomian's polynomials (see Seng, Abbaoui and Cherruault 



1996 , Abbaoui, Cherruault and Seng |1995| ). 

Definition 1. We say that the FD-method for the Cauchy problem pD con- 
verges (to the exact solution of problem ^) on [x , x + H ) , < H < +00, 
if there exists a real number h > such that for any grid u> ^3) with h < h, 
series Q) ; with the terms computed from converges uniformly (to 

the exact solution of problem ^) on [xo,xo + H) . 

In the present paper we consider a particular case of the Cauchy problem 



1|) when r = 1, Li (•) = (•) , more precisely 



d 

dx 



u(x) — N (x, u (x)) u(x) = 4> (x) 



(9) 



u (x ) = m , x G [x , +00) . 
The main result of the paper is presented by the following theorem: 
Theorem 2. Let the Cauchy problem IM) satisfies the following conditions: 



1. iV (x, u) = Y^ a i ( x ) u \ where a, (x) G C ([x , +00)) . And there exists 

i=0 

a sequence of real numbers Bi > 0, . . . i 



0,1, . . such that 



sup |oi (x)\ < B { , < B t G M 

xS[:ro,+oo) 
00 

and i/ie series E convergent Vw G R; 

i=0 



0,1,..., 



2. (f) (x) is a continuous and bounded function on [xq,+oo), with 



sup \(p(x)\ = k < +00; 
3. (AT (2, u) u)' u = N (x, u) + uN' u (x, u) < -a, < 0, Vx G [x , +00) , 

T/ien /or any initial condition uq G R £/ie solution u (x) of problem exists 
on [xq, +00) . The FD-method converges on [xq, +00) to £/ie exaci solution of 
problem |$|). Moreover, the following error estimations hold true 



sup 

xe[o,+oo) 



w (a; — u (x, 



< 



C 



(h/R) 



m+l 



sup 

ie[o,+oo) 



M ( X ) — U \ X) 



< 



[l + m) 1 + s 1 -h/R 



if h<R, (10) 



,1+e 



i=m+ i (j + 1) 

were e, R, C are positive real numbers that depend on the input data of the 
problem only. 



2. Justification of the FD-method algorithm for solving the Cauchy 
problem on the infinite interval. 

To begin with, let us introduce some useful notations. For each real- 
valued function u (x) defined on (a, 6) , a < b < +00 we denote 

IM a Ollo,oo,<o,&) = SU P \ u ( x )\> 

xd(a, b) 

and Vw (x) G C ((a,b)) we denote 



M^) II l lDO) (a, b) 

d 

dx 



max<j ||w(s)||o i00l<aj6> , 

0,oo, (a, b) 

For the future convenience we would like to define such notations 



u (x) 



\ U ( X )\\i,{a,b) = \\ U ( X )\\i,oo,{a,b) » \\ U ( X )\\i = \\ U ( X ) 



\i,oo, [x , +00) ' * 1' ^- 

(12) 

Before proceed to the proof of theorem [2] we need to prove several auxiliary 
statements, presented below. 
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Lemma 3. Let conditions 2) and 3) of theorem^ be fulfilled, and N (x, u) E 
([^Oj +oo) x R) . Then for any hi > 0, i — 1, 2, . . . i/ie solution (x) 



of the following Cauchy problem with piecewise constant argument (see. Akhmet 



d 

dx 



(x)—N (x, V, (x) — (j){x) , X E [Xi-i, Xi) , Xi = Xi-i + hi 



u {0) (x ) 



u E M, [u 



(0) 



x 



0, i = 1,2,... 



(13) 
(14) 



exzsis on [xq, +oo) and |m^ ^ (x) < fi = max {\uq\ , ^} , Vx E [xq, +oo) . 



Proof. Since equation (13) is linear on each segment Xj] and iV (x, u) 

is defined on the whole real line as a function of u, the solution of problem 
(13), (14) on the interval [xq, +oo) exists obviously: 



u 



(0) 



(x)=exp<^ / N{iM Q) {x^ l ))di\u^{x l ^) + 



Xi-l 



(15) 



+ / exp { f N (t, m (0) (xi_i)) dr f <j) (£) d£, x E [xj_i, X;] . 



Thus, it remains only to prove that this solution is bounded by some posi- 
tive constant /i. The conditions of this lemma together with the mean value 
theorem provide us with the following inequality 

N (x, u)u 2 =u [N (x, u)u — N (x, 0) • 0] = 

= u \{N' U (x, u) u + N (x, u)) u) < —au 2 , u E [0, u) . 

If we, additionally, take into account that the function N (x, u) , is continu- 
ous, we will obtain the inequality N (x, u) < —a, Vu E K. Then multiplying 
both sides of equation ( 13 ) (for i = 1) by (x) and using the above estimate 
for iV (x, u) we get the following 



1 d 

2dx 



( x )) 2 = N (x, u (0) (x )) (m (0) (x)) 2 + <P (x) u {0) (x) 



1 d 
2dx~ 



(u® (x)) 2 < -a («(°) (x)) 2 + k \u® (x)| . 



(16) 



Inequality (16) states that solution (x) of problem (13), (14) on [xo, Xi] 
can't leave a segment [— /i, /i] . 



6 



As a matter of fact, if we suppose that there exist G [x ,xi] such 

that x* < x* and 



a 



(x)\ > /i, Va; G (x*,x*\ 



(17) 



then we immediately obtain from (16) 

d 



dx 



(it<°> (z))' 



< 0. 



This contradicts our assumption (17). So, (x) G [—fi, /j] , Vx G [xo, a?i] . 
In a similar way we can prove that the solution of the Cauchy problem (13), 
(14) on [xi, xq\ belongs to the segment [— //, /i] too and so on. This completes 
the proof of the lemma. □ 

The following lemma can be considered as a generalization of lemma 2.1 
from |Gavrilyuk et al.| |2009 . 



Lemma 4. LetN(x,u) G ([x Q , +oo) x E) andu ij) (x) G C 1 ([x , +oo)) 
j — 0, 1, 2, Furthermore, let there exists a function N (u) G C 00 (R) siic/i 



sup 

a;S[a;o, +00) 



QP 



duP 



N(x,u) 



< 



dP 
duP 



N{\u\) Vu G R, \/p G N|J{0} 



T/ien t/ie following inequalities hold true 

\\A k (N(x,(-)); u^(x),...,u^ (X)) 



-A* (#(*,(■)); « (0) (^-i),---,« (fe) (^-i)) 



l0,[a;j_i, Xi] 



< 



< h,A k [N'(u)u- ||u {0) (x)|L r 



<h t A k (N'(u)u; Wu^ix)^ 



[xo, +00) ' ■ ■ ■ ' 



\ U lll,[a^_i, Xi) 

j I It ("^) || l^xo, +00) 



< 



Vx G [x^Xi] , fc = 0,1,2, 



.* = 1,2, 
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Proof. Throughout the proof we use the notation 



ftp 

iV(p) ( u ) = - — N (x, u) , u ij) = u ij) (x) , v®! = u ij) (Xi-l) ■ 



duP 



Let us fix any arbitrary I . It is well known that the follow- 

ing representation for the Adomian's polynomial holds (see |Abbaoui, Pujol, 
Cherruault, Himoun and Grimalt [2001 , Seng et al. [1996] ) : 



A k (N(-); «<°>,...,u«) 



(« (0) ) 



(ai-«a) r u (fc-l)j («h-i-«fc) [ M (fc)] a * 

(ai-a 2 )! "' (a fc _i-a fc )! (a*)! 



ai+...+a fc =fc 
ai>a 2 >...>a fe , aiSNUW, 

Therefore 



A k (N(.);uV\..., U M)-A k (iV(.); u2?i,-..,«2?i) 



< 



< 



r (i)i( a i-«2) r (fc)] 
V iv(-) (««») i^J _ . . . ^± 



0,[a;i_i, a*] 



Qi+...+afc=fc 



u 



(i) 
t-i 



(fc) 



i-1 



(ai-a 2 )! («fc)! 



"* E 

ai + ...+afc=fe 



■it 



(0) 



dwW ( 



dx 



Q,[xi-i, Xi] 

(1)1 ( a l- Q 2) 



[#*)] 



"A- 



x=5; (ai - a 2 )! ' ' ' (a*)! 



(ai — « 2 )! (a*. — 1)! dx 



< 



OJiEi-l, Xi] 



< 



E 



Qi+...+a fc =fc 



(ai - a 2 )! 



1/ 



(fc — 1) lic^fc— 1— <=*=fc |L ( (Jfc)||a* 

11 I I (Jj 11 



ll,[xi_l, Xi] 



(a fc _i-a fc )! («fc)! 



x 



x 



(0) )« (0) || nr ,+«! ||JV (0Il) (t* (0) )|| or 

/ M0,]xi_i,Xi] II V / 110, [Xi-i, Xi\ 



< 



8 



< 



a-i+...+a k =k 



^ 1 1 1 j — 1 5 3? j ] 

(«i - q 2 )! 



I 1 1 \ %i — 1 i 31 i 



I a*, 



X 



(a fc _! - a fc ] 



x 



)Ml u (o)| 



l,[xi_i, Xi] 



)||« (0) |L l + a l jV^)f||«(°)|| 1I .) 

J II lll,[xj_i, Xi\ ^ \\\ lll,[Xi_i, Xi\ J 



= h t A h (N> («)«; ||«<°) (x)|| lj[ei _ ii<Bi] , . . . , ||r*« WH^J < 

< hl A k (n> (u) u; \\u^ (x)\\ Uxo! +oo) , . . . , ||«W (x)|| Uxoj +oo) ) , 

where «™ = U U) (x) , 5? 6 [a^-i, Xj] , h* = x — And this is precisely the 
assertion of the lemma. □ 

The following lemma can be considered as a generalization of lemma 2.2 



from Gavrilyuk et al. 2009 



Lemma 5. For any scalar function N (u) G C°° ( 
holds true 

A j+1 (N{x); ^ o ,Vi,...,^,0 



i/ie following equality 



dP 



+i 



(j + 1)! 

oo 

where j = 0,1,..., f (z) = z%v i- 



I (N (f (z)) - (f (z) - V ) N' (V ) 



2 = 



i=0 



The proof of lemma [5] is trivial. 

Proof. ( Of theorem [^j Let the conditions of theorem [2] be fulfilled. It is easy 
to seen that for any initial condition uq G R, the solution of the Cauchy 
problem ^ exists and is unique on [0, +oo ) . Moreover , it be longs to the 

p.278). So, 



Demidivich 



1967 



segment [— /i, //],// = max{|M | , ^} (see. 

we are going to prove that this solution can be found by the FD-method. 
Let us fix any arbitrary infinite grid 

Q = {xi = Xi-i + hi, hi > 0, i = l,2,...; Xi — > +oo, i — > +00} , (18) 

h = sup hi 



9 



and embed the Cauchy problem ^ into the more general problem 
d 



dx 



u (x, t) — N (x, u t)) u (x, t) — 



— t {N (x, u (x, t)) — N (x,u (xi-i, £))} u (x, t) — (j) (x) , x £ [xj_i, Xi) , (19) 
[u (x, t)\ — 0, z = l,2,..., u (x , t) = u , Vt £ [0, 1] . 



Apparently, if we set t — 1 then the solution of (19) coincides with the 
solution of problem ^ : 

u (x, 1) = u (x) . 

If we set t = 0, on the other hand, we will get the base problem 
d_ u (o) ( x ) _ N (^ u (o) ( x ._^ u (o) ( x) = (x) , x £ [x 4 _ l5 Xi ] , i = 1, 2, . . . , 

«(°) (x ) = mo, [«<°> (x)] x=x . =0,i = l,2,..., 

(20) 

which is analogous to problem Q. Then, we assume that the solution of 
problem ( 19 ) can be found in the form of series 



u (x, t) = Y^ tju{j) i x ) > vt e [0, 1] , Vx £ [x , +oo) , 

j=0 



(21) 



d d 

— u{x,t) = J2 tj j- M (i) (x),Vt£ [0,l],Vx£ [xi_i, x^, 2 = 1,2,..., (22) 



where coefficients (x) depend on x only. By substituting (21) and (22) 
into ( 19 ) and comparing the coefficients in front of the powers of t, we obtain 
the following recurrence sequence of linear problems for (x) (with a 

piecewise constant coeficient): 



£ U V+V ( x ) - N (x, u<°) m^' +1 ) (x) = 

^ (x, «(°) (zi_i)) (x) m^' +1 ) (xi_0 + (x) 
x £ [xj_^, Xj] , i 1,2,..., 



(23) 
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where 



F (i+i) (x) = j2 A j+1 „ p (N (x, (■)) ; (xi-t) u«+^ (xi_i)) (x) + 
p=i 

+ t [Aj-p (N (x, (•)) ; u<°> (s) , t*« (*),..., (x)) - 
-A,-_ p (iV (x, (•)) ; (x^) , U (D , . . . , u V-p) (x) + 
(iV (x, (■)) ; u {0) (xt-i) , u« (xi-i) (xi-i) , 0) u(°> (x) , 
(x)] x=x _ = 0, * = 1, 2, . . . , u^' +1 ) (x ) = 0, j = 0,1,2, 



It is easy to see, that problems (23) are analogous to problems ([5]). Let us 
express equation (23) in the equivalent form 



<(x,«(°) (x,_0) 



-0 / |« (0) (0^+ 



(24) 



u^ +1 ) (x 



0, i = l,2,..., (x )=0, 



where 



Ql (x) = iV (x, u(°> (x^)) + < (x, (x 4 _!)) u<°) (x^O . 



The continuous solution of the Cauchy problem (24) can be represented in 



11 



the following form 



(x) 



exp \ j Qi (0 d£ } + 



+ / exp \j qi (r)dr \K(tuW(x i ^ 1 ))x 



Xi-l 



Xi-l 
X 



1 \ ^ J — l 



+ / exp <j / 9i (r) dr\ (£) d£, x G [x^, Xi ] . 



Xi-l 



(25) 



On the other hand, from (23) we will obtain 



u J ' (x) 



exp { / m (0 d£ } + 



+ J exp{jn t (rj) drj N' u (£, ti<°> (x^)) u<°) (0 



It 



(i+i) 



>i_i)+ (26) 



+ / exp{jn t (v)dv\F^(Odt, 



Xi-l 



where 



rti (x) = N (x, v, (xj_i)) x G Xj] . 

If we differentiate (|26|) with respect to x we will obtain 



12 



dx 



U yj ' [X] 



m (x) exp | / m (0 d^j + ^ (z, u<°> (Zi_i)) u(°) (x) + 
+ ni (x) J expi/ni^drliVifcuW^i-i))^ )^)^ 



(x) / exp <( / Tij (r) dr \ F^+V (f ) ^ + ( 



(27) 



Equation (27) could be expressed in the another form 



dx 



+m (x) J exp I J m (r) dr \ (f) d£ + (x) 



(28) 



where 



Pi (x) = ni (x) exp I j m (0 dO + < (a?, u {0) (x;-i)) u<® (x) + 

a; ( i I 

+rn (x) J exp I J m (r) dr \ (f , (ar^j)) w (0) (£) d£, x G a? 4 ] . 

<m~i U J 

(29) 



From (28), (29) and lemma |3| we obtain the following estimates 

"- Uij+1) ^klx^x, < lift WUo,^,*] M\ + 



I da: 



+ (1 + Nhe Nh ) (x)|L r iV = max |iV(x,u)| 

V 'II v ' \\0,[xi-i, Xi] |„|< M 1 

a:S[a;o,+oo) 



(30) 



for i = 1,2,.... For abbreviation, we denote (using the result of lemma |3 

5= max (sc, m)| (\N (x, u)\/i + k) < +oo, 

\u\<n 
x£[xq.+oo) 

C= max |iVy (x, u)\ // < +oo. 

\u\<fi 
x£[xo,+oo) 
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Then, from (25) we obtain the estimates 



lk (i+1) (*)IL i< 

Chi d 



1 + 



\u 



(x<-i)| + 



—vV +r > (x) 
ax 

< Xj )\ < 



+ hi\\FV +x > 



x 



Chi 



dx 



u 



u {j+1) (x) 



0,[xi-i, Xi] 



0,[x , +00) 



Bh3 



u 



+ h l \\F (j+1) ( 



i — 1 2 
lo,[x , +00) ' ' ' 



Combining the last two inequalities, the result of lemma 3 and (30), we get 



I^MIo,^,^] ^ [ l + Bh *] h 0+1) (^-i)|+^||^ (i+1) (^)||o. 



u 



were 



B = ^ + ^,D = l + h~(l + Nhe Nh ), 



[xo, +00) ' 

(31) 



0,[xo, +00) ' 

(32) 



p = Ne Nh + C + Nhe Nh C = C + Ne Nh (1 + hC) > \\p (x) \\ 0t[xo>+oo) ■ 
Without loss of generality, we can assume that (^) llo,[a:o,+oo) 

> 0. 

Since that we denote 



IO,[a;o, +00) 



and rewrite ( 32 ) in the following form 



y . < [ e -a^ + Bti:] + hiD, y = 0,i = l, 2, 
It is easy to see that for all hi > 



(33) 



-ahi 



+ Bti <l-ahi + ti [B + 



a 



l-hi[a-hi[B + 



a 



Thus the inequality 



hi < h< 



a 



2B + a 2 
14 



z = l,2,..., 



(34) 



implies 



- ahi + Bh? < 1 "' 



hiOL 

2~" 



The last inequality, under conditions (34), guarantee that the recurrence 
sequence 

(35) 



Y, = ( 1 h/^ ) ' \ YD. i 1.2 V„ = 0. 



is a majorant for the sequence yi (33). 

It is easy to seen that (35) can be reduced to the form 



where 



Zi = (l - hi) i = 1,2, Z = 1, 



h — h — 7 — 1 _ V V — _1 7 — 19 

<H 'H r\ ' i 1 i 1 i o'n ' ' ' ' ' 



hence that 

i 

Z i = Y[(l-h p ), i = 1, 2 Z = l. 

p=i 

We conclude from the last expression that the inequalities 



< hi < fii = min 
implies the estimates 



ft 25 + ft 2 



Vi e N 



(36) 



From now on we orient the grid cD (18) by requirement that (36) is satisfied. 
Thus the inequalities (30) and (31) yield 



\u {j+1) (x) 



\0,[xi 



X 



l0,[xo,+oo) 



dx 



u u+1) (x) 



< (pD- + — + 1 ) exp 
\ ft ft ' 



|_pO'+i) 



Combining the last two inequalities, we obtain the estimate 

II v ' \\l,[x ,+oo) — II v / M0,[ao,+oo) ' 



\0,[x ,+oo) ' 



(37) 
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where 



a = max j ^(l + Bfij) —D + HiE^j , (pD- + — + exp 



To make the following estimations we need to use the results of lemmas 
4] and [5] The requirements of these lemmas are satisfied if we set (in the 
notations of lemma [4]) 



n(u) = ik i^Wo^' <J2 BiVji ' Vn G R - 

i=0 i=0 



(38) 



(37) is equivalent to 



O+i) 



x I, < 



<<MX>+i-* (#(«); \\u^ (x)^ , . . . , 



UO+i-p) f x ) 



\u ip) (x)\l + 



p=i 
j 



+h^2Aj- p (N' (u)u; ||m (0) (x)^ , . . . , \\u {j - p) (x)|| x ) ||« (p) (x)\\ 1 + 

p=0 



+ 



(x 



O' + l)! 



HE 



z s lk (s) (x) 



,s=0 



s=l 



j = 0,1,2,... 



z=0- 



Denote 



J = 0, 1, 



3 ~~ W 

Let us define a sequence {V^}^ by the following recursive formula 



u) ; Vn, . . • , K- 



, V , ■ • • , Vj+l-p j Vp 



v v + 



,p=l 



(39) 



(40) 



(41) 



+ Aj. p (n' {u)u;V ,..., Vj-pj V p + 



p=0 
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V j+1 V N' (V ) ,j = 0,1,.... 



,s=Q 

or, in more convenient form, 



2 = 



(7 



1 + aV N' (Vb) 



£ (iV (it) ; Vq, . . . , ^+i_p) r p + (42) 



. p=0 



+ ^ (iV' (u) u;V ,..., Vj- P ) V p f . 

p=0 

where Vq — fM. 

It is easy to seen that v i < Vi, Wi € N |J {0} . 

We are now in a position to prove that for h (18) sufficiently small the 
assumptions (21) and (22) hold. In order to do that, it is sufficient to show 

(43) 



that the function 



3=0 



has a nonempty open domain. 
From 

g(z)-v 



From (42) we obtain 
a 



1 + aV N> (V Q ) 



N(g(z))-N(V ) + zg 2 (z)N'(g(z)) 



(44) 



Let us express z from (44) as a function of g : 

1 



9 2 N'(g) I s 
V < 9, S 



(<?-V )- [N(g)-N(V ))g 



a 



1 + aVoiV' (V 



(45) 



Evidently, the function z (g) (45) is defined and continuously differentiable 
in some open neighbourhood of the point g = Vq. To prove the existence of 
the inverse function g = g (z) defined in some open neighbourhood of the 
point z = it is sufficient to show that z' (Vo) > 0. The last fact immediately 
follows from formula (45): 



z' (Vo) = lim 



z{g)-z(V ) 



s->vb g — Vo 
17 



lim 



smv V 2 N'(V ) V E 



aV 2 N> (y ; 



N(g)-N(V ) 

g-v ; 
> o. 



(46) 



From (38) it follows that N (u) — > +oo as u — > +oo. Hence, from equation 
(45) we have: 

lim z (g) < 0. 

g— >+oo 



The last inequality together with (46) provides the existence of a point g max 
such that g max > V , z' (g ma x) = and \/g e (V , gmax) Z 

(g) >o. R = 

Zmax Z {gmax) is the radius of convergence of the power series M3b. In 
other words, 



R J Vj < C- 



1 



(j + l) 1+£ ' 

for sufficiently small e > 0, where C is some constant independent of hi, i 6 
N. 



Returning to notation (40), we obtain 

C 



\u (j) (x) 



I l,[a;o, +oo) 



< 



U + 1) 



1+e 



,J = 0,1, 



(47) 



This estimate gives us the following sufficient condition for series (43) to be 
convergent: 



h 
R 



< 1. 



Finally, we have proved that for h > sufficiently small (to be exact, 
h < min{/ii, R}) assumptions (21) and (22) are fulfilled. Therefore, it is 
remain to prove that the sum of the uniformly convergent series rt2J) is a 
solution to problem rt9|). To do that we need to add up the base problem (20) 
and equations (23) for j = 0, 1, .... As a result we obtain 

oo , oo j 

E t uU) E V* ( x ' (•)) ; u(0) (*),■■•, «°'" p) (*)) « (p) (*) = * (*) 



j=0 



j=0 p=0 



(48) 
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x G U Xj) . Using the theorem about the composition of two power 

i=i 
series (see 



Fihtenholts 1966 , p. 485), it is easy to prove the equality 



.V ( x^t l u {i) (x) J2 fu(l) ( x ) 

i=0 / i=Q 



= E tJ E (*> (•)) > m(0) (*) > • • • > (*)) (*) . 

i=0 p=0 

Vt G [0, 1] , Vi G [xq, +oo) . Thus, taking into account the uniformly con- 



vergence of the series Yl ~h u ^ ( x ) on eac h interval (xi-i,Xi 

i=o x 



1,2, 



we can rewrite ( 48 ) in the following form 
d 



dx 



u(x) — N (u (x)) u(x) = <p (x) , Wx G (xi-\, Xi) , u (x) — 

i=o 



i=l 



(49) 



Recall that u (x ) = u . Then the continuity of u (x) , the existence and 
uniqueness of the solution u (x) of Cauchy problem ^ on [x , +oo) along 
with (49) imply our final goal u(x) = u (x) , Wx G [xq, +oo) . This means 
that the FD-method for the Cauchy problem ^ converges to the exact so- 
lution of the problem in the sense of definition [T] 



We leave it to the reader to verify the error estimates (10), (11). They 

□ 



could be easily derived from (47). 



3. Examples. 

Example 1. As an example let us consider the following Cauchy problem 
d 



dx 



u (x) = —u (x) — u (x) + cos (x) + sin (x) + sin (x) , u (0) = 0. (50) 



It easy to seen that the exact solution of problem (50) is u* (x) = s'm(x 



To solve problem (50) numerically, we, first of all, may try to apply the 



ADM. We approximate the solution by the m— th partial sum ua of the series 



(x) , where (x) can be found from the sequence of Cauchy problems 

i=0 



.0 



d 

dx 



uf (x) 



-u 



(0) 



x) + cos (x) + sin (x) + sin 3 (x) , (0) = 



19 



.(0 



d 



dx 



u A (x) = --t 



(x) 



i=0 



^(z), < ; (0) = 0, z = l,2, 



.(0 



t=o 



The results are presented on the Fig. 1. The graphs on Fig. 1 show that the 



if 

-5 
-10 
-15- 
-20 



U(X) 




a) 



b) 



Figure 1: Example 1. ADM, continuous line: u* (x) — sin (x) ; o : ua (x) ] D : U A {%) \ 

2 3 

: u A (x) ; + : (a;) ; 



ADM for the Cauchy problem (50) is divergent on [0, 6] 



The application of FD- method to problem (50) with the grid 

1 



Xq — 0, Xi 



1 . . . , 144 



yields results which are presented on Fig. 2, Fig.3. Here we use the notations: 

5i (x) = sin (x) — u (x) , i — 0, 1, 

The graphs on Fig. 2 and Fig. 3 confirm the exponential convergence rate 
of series (j2j) to the exact solution of problem (50). 
Example 2. 

Let us consider the Cauchy problem 



d u (x) + ( -4= + 1J u 3 (x) = ^-^= + 1J sin (2y^ + x) 



dx 



x 



(51) 
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a) 



b) 



Figure 2: Example 1. FD-method. The graphs of absolute errors, a) - 5q (x) ; b) 
Si (x) . 




0.0005 - 



.' 
0.0005 : 
-0.001 : 



Y 50 



a) b) 

Figure 3: Example 1. FD-method. The graphs of absolute errors, a) 
S 3 {x);xe [0, 48] . 



5 2 {x); b) 



x e [0, l] , u(o) = l. 



Problem (51) has a singularity at the point x = : lim -j-u(x) = +oo. 

i->0 
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Figure 4: Example 2. FD-method. + : u (x) , □ : u (x) , o : u (x) . 




a) b) 
Figure 5: Example 2. FD-method. Graphs of discrepancy. a) - v a (x) , b) - V\ (x) . 



Using routines from the computer algebra system Maple 12, we tried to find 
the solution of problem (f5Tj) either analytically or numerically, but all that 
was in vain. In the present version Maple is unable to solve such problems. 
The conditions of theorem [2] are not fulfilled for this problem, either. Re- 
gardless of that, the ideas of ADM and FD-method are naturally applicable 
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Figure 6: Example 2. FD-mcthod. Graphs of discrepancy, vi (x) . 



to this problem. However, as it turned out, the ADM is divergent on [0, 1] 
in this case. 

To apply the FD-method we need to introduce a grid on the segment 
[0, 1] 

u = {0 = x , Xi = 0.05i, i = 1,2, ...,20}. 
The base problem is stated as follows 



d 

dx 



«(°> (x 



+ ^-±= + lj ( u (°) (xi-i)) 2 ^ (x) = (l= + lj sin (2^ + x) , 



x e m (0) = l. 

It admits the analytical solution 



(o) 



(x) = exp - 2^1 + (^ (0) « (0) fo-i) + 



+ 



| + l) exp ( (2v/f + r) (« (0) (Xi-i)) 2 ) sin ( 2> /l + d£ 



exp (2^ + ^) 



? * («W(x i -i)) 2 )« (0) (x i _ 1 ) + 
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Similar analytical formulas were obtained for (x) and (x) . For error 
control we use the discrepancy 




The results are presented on Fig. 4 - Fig. 6. As above we use the notation 
u ( x ) = Yl u ^ ( x ) ■ ^ i s eas y to seen that the discrepancy is decreasing 

i=0 

exponentially even behind the point of singularity x — 0. 
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